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Spatially inhomogeneous electron distribution around a 
single vortex is discussed on the basis of the Bogoliubov-de 
Gennes theory. The spatial structure and temperature de- 
pendence of the electron density around the vortex are pre- 
sented. A relation between the vortex core charge and the vor- 
tex bound states (or the Caroli-de Gennes-Matricon states) is 
pointed out. Using the scanning tunneling microscope, infor- 
mation on the vortex core charge can be extracted through 
this relation. 

PACS number(s): 74.60.Ec, 74.25. Jb, 74.72.-h, 61.16.Ch 

Electric charging phenomena around vortices have the 
potential of becoming one of the key features in the 
physics of the mixed state in type-11 superconductors. 
Until quite recently, little was known about the electric 
charging inherent in vortices, while it has been well recog- 
nized since the 1950's that each vortex line carries a quan- 
tized magnetic flux. Only recently, it was noticed that an 
electric charge accumulates around a static vortex line in 
type-II superconductors. Khomskii and Freimuth,EJ and 
Blatter et alu theoretically discussed the electric charg- 
ing around a vortex.cl If the electric charging of vortices 
is experimentally confirmed, it will open the door to an 
unexplored field in which one expects various electromag- 
netic phenomena to originate from the electric charge 
trapped by vortices. 

In spite of the growing interest in vortex core charg- 
ing, firm experimental evidence of the charging is lacking 
at the present. However, various experimental attempts 
are now in progress and are on the verge of detecting a 
charge accumulation inside vortex cores. One such ex- 
periment is a polarized neutron scattering investigation 
of the flux line lattice in Nb by Neumann et ala They de- 
tected a nonzero nuclear contribution to the Bragg peaks 
corresponding to the periodicity of the flux line lattice. 
This experimental result strongly suggests existence of 
the vortex core charge; if each vortex constituting the 
flux line lattice traps electrons or holes, the response of 
the underlying nuclei to these will indue©, a distortion in 
the nuclear lattice around the vortices au Various types 
of experiments which will attempt to detect the vortex 
core charge are also planned. In addition, an experiment 
to observe the temperature T and magnetic field H de- 
pendence of the vortex core charge is expected in order 
to establish the existence of the vortex core charging. 
Therefore, it is certainly desired that detailed theoretical 
predictions for the temperature or magnetic field depen- 
dence of the vortex core charge should be presented for 
experimental verification. 



In this paper, we present the structure of the carrier 
density around a static single vortex and its temperature 
dependence, solving self-consistently the Bogoliubov-de 
Gennes (BdG) equation. On the basis of the solution of 
the BdG equation, we discuss not only the temperature 
dependence but also the relation between the charging 
of the vortex core and the so-called Caroli-de Gennes- 
Matricon (CdGM) states (or the vortex bound states). 

The CdGM states, i.e., low-energy excited states due 
to vojjticesja were first discussed theoretically by Caroli 
et alu Their existence was experimentally confirmed by 
Hess et al.p who observed spatial dependence of the ex- 
citation spectra around a vortex with scanning tunneling 
microscopy (STM). The local density of states (LDOS) 
around a vortex, probed by STM, depends on the Bo- 
goliubov wave functions of the CdGM states Uj(r) and 
Vj(r), labeled by the quantum number j. The LDOS 
N(r,E) (to be exact, thermally smeared LDOS, i.e., the 
tunneling conductance) is given as 

N(r, E) = Y1 [Mr)| 2 /'(£ - E 3 ) + \v 3 {v)\ 2 f'{E + Ej)] , 
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(1) 

where Ej is the eigenenergy and f(E) the Fermi function 
(the prime represents the derivative). The STM enables 
us to extract detailed information on the wave functions 
around a vortex. Here, we notice that the carrier den- 
sity around a vortex, n(r), also depends on these wave 
functions: 

n(r) = 2 J2 [K(r)| 2 /(^) + h(r)| 2 {l-/(^-)}]. (2) 

Ej>0 

The electric charging (or the inhomogeneous electron 
density distribution) around a vortex is related to the 
LDOS through the wave functions Uj(r) and Vj(r). This 
suggests unique potential ability of the STM; the struc- 
ture of the LDOS probed by STM relates to the spatial 
structure of the vortex core charge. 

Regarding the previous theories of the mechanism of 
the vortex core charging, Khomskii and Freimuthlil based 
their scenario on a normal-core model. Assuming that 
the vortex core is a region of normal metal surrounded 
by a superconducting material, they considered that the 
corresponding difference in the chermcal potentialcUeads 
to a redistribution of the electrons. El Blatter et aln dis- 
cussed the charging mechanism, considering spatial vari- 
ation of the pair potential A(r) around a vortex. On 
the basis of the zero-temperature version of eq. (|^) , they 
obtained n(r) by combining the spatial variation of the 



1 



wave function v(r) with particle-hole asymmetry in the 
normal-state density of states at the Fermi level. The 
discussion was, however, based on a wave function which 
had the same fomi as the uniform solution of the BdG 
equation, namelycl 



N 



E k = 



|A(r)P 



(3) 



The spatial variation of Vk (r) was directly determined by 
the local value of A(r), which is not exactly appropriate 
for the vortex system. It is desired that one should base 
the calculation on the exact wave functions of the CdGM 
states. 

Prompted by this motivation, we will self-consistently 
solve the BdG equation to obtain the exact wave func- 
tions Uj (r) and Vj (r) of the CdGM states (including the 
extended states above the gap). We start with the BdG 
equation given, in a dimensionless form, by 
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li Uj(r) + A(r)«y(r) = E jU] (r), 



v j (r) + A*(r)u j (r)=E j v j (r), (4) 



where fj, is the chemical potential and £o( =ii f/Ao) is the 
coherence length [A is the uniform gap at T = 0, and 
/cf (vf) is the Fermi wave number (velocity)]. In eq. (Q), 
the length (energy) scale is measured by £n (Ao). For an 
isolated single vortex in an extreme type-II superconduc- 
tor, we may neglect the vector potential in eq. (Q). To 
maintain macroscopic charge neutrality in the material, 
in eq. (|^) we constrain the electron density in a uniform 
system to be constant on the temperature. We use /i de- 
termined at each temperature by this constraint, which 
is equivalent at zero temperature to eq. (4) of ref. [| The 
pair potential is determined self-consistently by 



A(r)=g £ Uj (r)v*(r){l-2f(Ej)}, 



(5) 



where g is the coupling constant and lod the energy 
cutoff, which are related by the BCS relation via the 
transition temperature T c and the gap A . We set 
uju = 20 A . We consider, for clarity, an isolated vor- 
tex under the following conditions, (a) The system is 
a cylinder with a radius R. (b) The Fermi surface is 
cylindrical, (c) The pairing has isotropic s-wave symme- 
try. Thus the system has cylindrical symmetry. We write 
the eigenfunctions as uj(r) = u Tli ;(r) exp [i(l — \)0\ and 
Vj(r) = v n j(r)exp[i(l + ±)6] with A (r) = A(r)exp[-i0] 
in polar coordinates, where n is the radial quantum num- 
ber and the angular momentum |/| = i ~, |,- ••. We ex- 
pand thej-eigenfunctions in terms of the Bessel functionsS 
J m {r) asEl 



1-41(0. 



v„,i{r) =^2<i ni (p ill+ x l (r), 



(6) 




where 4> im (r) = [V2/RJ m +i( )]J m (a im r/R) and 
is the i-th zero of J m (r). We set R = 20£ - The BdG 
equation is reduced to a 2N x 2N matrix eigenvalue prob- 
lem. This useful technique to solve eq. (§), develop^ 
Gygi and Schliiter,E3 has been utilized in some cas 
Our system is characterized by a parameter kp^Q^t^i im- 
portant for the present problem. From our standpoint, 
all interactions between the quasiparticles are renormal- 
ized to g in eq. (|^) and additional screening does not 
exist in the Hamiltonian. The screening for the charge 
ordering is excluded as in the charge density wave stud- 
ies.lia If some screening effect is considered, in principle 
we may take it into account as an external potential in 
eq. ([|) and solve self-consistently the equations together 
with an additional equation, e.g., the Poisson's equation. 
Such a study, if meaningful, is left for a future work. Us- 
ing the calculated Uj(r) and Vj(r), we obtain the LDOS 
iV(r, E) and the carrier density n(r) from eqs. ([!]) and 
(||), respectively. 

In Fig. Fy, we present the spatial structure of the carrier 
density n(r) around the vortex at several temperatures. 
The Friedel oscillation appears at low temperatures, be- 
cause each wave function of the low-energy CdGM states 
oscillates with a period ~ fcp 1 . It is striking that the 
carrier density at the vortex center exhibits strong tem- 
perature dependence and leads to a substantial charging 
at low temperatures. 




= 0.0 
0.3 
0.6 
0.9 

(extended states) 0.0 



FIG. 1. The spatial variation of the carrier density n(r) (in 
arbitrary units) around the vortex for several temperatures 
and fcp^o = 4. The contribution of the extended states to 
n(r) at T = is shown by the dotted line. The difference 
between this contribution and the total n(r) at T — gives 
the contribution of the bound states. 

The carrier density at the vortex center in Fig. |l| de- 
creases with respect to that far from the core. Conse- 
quently, in the case of the present electron system (i.e., 
the two-dimensional free electron system), the sign of the 
vortex core charge is opposite to the sign of the electron 
which is the dominant charge carrier in the present case. 
When the dominant charge carriers are holes, we only 



2 



have to treat these holes as carriers in that system in- 
stead of the electrons and there are no changes in the 
formulation [eqs. (@), (§), (@)-(@)]. The density of the 
dominant carriers (holes) decreases near the vortex cen- 
ter in this case as well. 

The density of the dominant carriers decreases near 
the vortex center, as long as the wave functions arouiad. 
a vortex for the dominant carriers are given by eq. (g) .E_ZJ 
This is related to particle-hole asymmetry in the LDOS 
inside the vortex core and can be understood in connec- 
tion with the CdGM states as follows. In the defini- 
tion of the angular momentum I in cq. (^]), the bound- 
state energy spectrum is E\ > for I > (and Ei < 
for I < 0), where E\ = —E-i. In Fig. z, we show 
the spectral evolution obtained from eq. (1). In sys- 
tems where &f£o is small (the quantum limit) , the asym- 
metry in the LDOS appears conspicuously. The two 
largest peaks near E = are noticeable [the peaks A 
and B]. The peak A at E = E l=i / 2 (> 0) is composed of 
\ui=i/2(r)\ 2 (= h=-i/ 2 M| 2 ). The peak Bat£ = E_ 1/2 
(< 0) is composed of |t;jy 2 (r)| 2 (= |tj,_ 1 / 2 ( r ) | 2 ) • From 
eq. (||), ^1/2(0) ^ and vi/ 2 {0) = because J m (0) 
only for m = 0. The asymmetry between Ui/ 2 (r) and 
Vi/ 2 (r) leads tojjae particle-hole asymmetry in the LDOS 
inside the core.E3 Now, according to eq. (||), n(r) is con- 
structed from the wave functions which belong to E > 0. 
The contribution from the extended states (E > Aq) is 
presented as the dotted line in Fig. E[ The remaining 
contribution to n(r) comes from the bound states. The 
lowest bound state v 1 / 2 (r), which belongs to the lowest 
bound state eigenenergy E\j 2 > 0, predominantly deter- 
mines the structure of n(r) in the vicinity of the vortex 
center. The amplitude |wi/2(^)| 2 is equal to that of the 
peak B in the LDOS. The spatial profile of n(r) is deter- 
mined by the shape of |vi/2(?")| 2 , i-e., the peak B. Since 
l w i/2( r )| 2 decreases to zero with r — > as seen from the 
spatial profile of the peak B in Fig. H, we can infer that 
ii{r) decreases near the vortex center. 
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FIG. 2. The spectral evolution of the LDOS N(E, r) (in 
arbitrary units) at T/T c = 0.05 and kp^o — 4. 

According to discussions! based on eq. ([|), the carrier 
density near the vortex center has a sensitive dependence 
on the slope in the density of states. It might be ex- 



pected that if the derivative of the density of states is 
negative, the carrier density increases at the vortex cen- 
ter. To examine it, we have investigated the case of the 
energy band, k 2 /2m + fc 4 /4m 2 £o, (see ref. which has 
a negative derivative of the density of states in two di- 
mensions. In the calculation with a fixed fi, the density 
far from the core certainly decreases with the growth of 
the gap A(T) on lowering T, which is consistent with the 
precondition of ref. ^. In this situation, on the basis of 
eq. (||), the density n(r) is naively expected to recover to 
the normal-state value on approaching the center r = 
where A(r) — 0. n(r) is then expected to increase at 
the center. However, according to results of the calcu- 
lation based on the wave functions of the CdGM states, 
n(r) decreases at the vortex center. We conclude that, 
the carrier density near the vortex center is determined 
by the electronic structure inside the vortex core, which 
is insensitive to the slope in the normal-state density of 
states at the Fermi level. 
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FIG. 3. The temperature dependence of the carrier den- 
sity no = n(r — 0) at the vortex center. In the figure Sno/rioo 
are plotted for several Aif^o, where 5no = \no — noo| and rioo 
(= rii n f) is the plateau density far from the core. 

Let us focus on the magnitude of the core charge. The 
carrier density at the vortex center, from which the or- 
der of magnitude of the core charge is estimated, exhibits 
substantial temperature dependence as shown in Fig. 0. 
We plot (5no/«oo, where 5no = \rio — n.oo|, no = n(r = 0), 
and Tioo is the plateau density Tioo — ?7-(-R/2), to which 
the calculated n(r) settles away from the core. We note 
that the &f£o dependence of the density, Sno/noo ~ 
(/cfCo) - " — (Ao/ef)", varies with the temperature (ep is 
the Fermi energy). Our numerical data show that a ~ 1 
near T = and a sa 2 near T = 0.5T C . The exponent a 
is crucial to the magnitude of the core charge. In most 
conventional superconductors, the parameter fcp£o is of 
the order of 100. It can be 1 - 10 in high-T c cuprates. 
Depending on the estimate of a, there can appear sub- 
stantial difference in the evaluation of the magnitude of 
the core charge. According to our results, a depends on 
the temperature as above. To estimate the total core 
charge Q v per unit length along the vortex axis, we con- 
sider the charging volume in Fig. [j] to be a cone with a 



3 



height Stiq and a base radius r% (2k-pri = n). n(r) almost 
recovers to initially at r% ~ ftp 1 at low temperatures. 
Q v is evaluated as Q v f» e7rrf<5no/3. We consider a pan- 
cake vortex in a layer, and the distance between each 
layer is d. In this case n^ = 27rfcp(27r/d)/87r 3 . We then 
obtain Q v <~ e(fci?^o) a ^ _1 at low temperatures. 

We should comment on the vortex dynamics in the 
context of the above temperature dependence of no, al- 
though the issue concerning the adynamics is seriously 
controversial at the present timeO Feigel'man et alxQ 
proposed a nondissipative transverse force acting on a 
vortex originating from Sno (see also ref. |l8| ). Kopnin et 
aLr-,reported that the effect proposed by Feigel'man et 
alx2 can be understood from the viewpoint of the spec- 
tral flow theory whep, no is regarded as the spectral 
flow parameter Co.E3'E3 The parameter Co is indepen- 
dent of the temperature. Hence it appears to be incon- 
sistent with the temperature dependence of no presented 
in this paper. Even in a neutral system with a fixed /i, 
no exhibits substantial tenmerature dependence in our 
calculations. While Kopniro discussed the temperature 
dependence of that force, the temperature dependence of 
6no itself at the vortex center seems not to be explicitly 
included there. We hopft,for a further investigation based 
on the CdGM solutionsQ to reveal possible mutual rela- 
tions between these theories (refs. £(J^l] and ^) and the 
significant temperature dependence of n(r) in the present 
paperO 

We point out a relation between the present work and 
STM -experiments. Maggio-Aprile et al£B and Renner 
et alEB observed spectral evolutions of the LDOS in- 
side the vortex cores in the high-T c cuprates. They de- 
tected particle-hole asymmetry in the LDOS near the 
core center (see Fig. 2 in ref. £6|). We expect that the 
asymmetry observed in the experiments has the same 
origin as the asymmetry shown in our Fig. || (see also 
ref. |l^) . We speculate that even if the superconductivity 
in the compounds consists of the preformed pairs or is in 
the crossover region between the BCS superconductiv- 
ity and the Bose-Einstein condensation, the Bogoliubov 
wave functions would still be defined. If so, the electronic 
state of the vortex core in_ihe compounds is understood 
as the Andreev scatteringE3 and it is the coherent state. 
From our results based on the Bogoliubov wave functions, 
we conclude that the particle-hole asymmetry inside the 
vortex core observed in the experimentsE3n3 implies the 
corresponding existence of the vortex charging. Accord- 
ing to another STM experiment by Renner et al.jEB the 
coherent electronic structure inside the core, observed as 
sharp structure of the LDOS, is smeared gradually by 
impurity doping. We predict that the vortex core charge 
decreases by impurity doping, because the charging is re- 
lated to the sharp LDOS structure inside the vortex core 
in our scenario. 

In summary, we investigated the electron density 
around a single vortex on the basis of the BdG theory. Its 
temperature dependence was presented. We expect that 
experimental data regarded as the vortex core charge will 



exhibit the temperature dependence as shown in Fig. |3|. 
If such dependence is observed, those experimental data 
will become solid evidence of the vortex core charging. 
We discussed the microscopic charging mechanism, which 
is independent of the slope in the normal-state density of 
states at the Fermi level, by considering the CdGM states 
around the vortex. We pointed out the relation between 
the vortex bound states, probed potentially by STM, and 
the vortex core charging, based on the inherent particle- 
hole asymmetry inside the vortex core originating from 
the CdGM states of the vortex. 
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